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2^ ' Abstract 

^ I A deformed differential calculus is developed based on an associative ^-product. In two 

Q^' dimensions the Hamiltonian vector fields model the algebra of pseudo-differential opera- 

■ tor, as used in the theory of integrable systems. Thus one obtains a geometric description 

i of the operators. A dual theory is also possible, based on a deformation of differential 

CN ] forms. This calculus is applied to a number of multidimensional integrable systems, 
such as the KP hierarchy, thus obtaining a geometrical description of these systems. 

^ . The limit in which the deformation disappears corresponds to taking the dispersionless 

CN ' limit in these hierarchies. 

^ I To appear. Journal of Geometry and Physics. 



Oh! 



1. Introduction 



Is there a common structure behind all integrable systems? There are many different types 
of integrable systems; (1 + 1) and (2 + l)-dimensional evolution equations (such as the KdV 
^ I and KP equations), chiral and harmonic map equations, integrable dynamical systems 
d ' (such as the Halphen and Kovalevskaya top equations), integrable nonlinear ordinary 
differential equations (such as the Chazy and Painleve equations), for example, and all 
these have various properties associated with their integrability (see, for example, [AC]). 
However, there is very little in the way of general theory, where the apparently disparate 
properties of various individual integrable systems could be understood in a consistent 
and coherent way. Indeed, there is not even a universal definition of what integrability 
actually is. 

One idea, proposed by Ward [W], is that such system may all be obtained from the 
(anti)-self-dual Yang-Mills equations and their generalisation by a process of dimensional 
reduction. For example, the KdV, NLS, sine-Gordon and Liouville equations may all be 
obtains from the (anti)-self-dual Yang-Mills equations with SL{2, C)-gauge groups, the 
only difference being the choice of symmetry group and spacetime signature [W,MS]. The 
key idea is not so much the self-dual Yang-Mills equations themselves, but the existence 
of a Penrose transform for such fields. Under such a transform the fields 'disappear' into 
the holomorphic geometry (certain holomorphic vector bundles over regions of an auxil- 
iary complex manifold known as twistor space) [W77]. More prosaically, this provides a 
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geometric framework for the Riemann-Hilbert problems used in the construction of solu- 
tions to the duality equations.The existence of this transform has been conjectured to be 
behind the idea of integrability. 

The paradigm (and the original example of such a transform) came from the (anti)-self- 
dual Vacuum equations [Pen]. The following formalism is due to Gindikin [G]. Consider 
the following system of first order equations depending on a parameter r = {tq, ti} G : 

where a;] are 1-forms. Let Q'^{t) be the bundle of 2- forms 

0'=(t) = u;\r) A u;\r) + ... + u;^P-\r) A u;^P{r) 
satisfying the conditions: 

• the {p + 1)*'' exterior power of vanishes; 

• the p*^ exterior power of Cl'^ is non-degenerate; 

• dn''(T) = . 

The bundle of forms then encodes the integrability of the original system. In the special 
case k — 1 ,p — 1 the metric defined by 

has vanishing Ricci tensor and (anti)-self-dual Weyl tensor. The form fl is related to 
various structures on the corresponding curved twistor space. 

One major problem with this geometric approach to the understanding of integrable 
systems is to find how systems such as the KP equation 

{4:Ut — 12uUx — Uxxx)x — ^'^yy 

fit into the scheme. There have been various attempts, some erroneous, to do this, the 
problem stemming from how to give a geometrical description to the pseudo-differential 
operators used in the derivation of the KP equation and its hierarchy. However, though 
there are problems with the KP equation itself, these problem vanish in a particular limit 
(the dispersionless limit) of the KP equation and one obtains a geometric description of 
this limiting case. Explicitly, let 



X = ex , 

Y = ey, 

T = et, 

U{X,Y,T) = u{x,y,t), 



then the KP equation becomes, in the hmit e — > the dispersionless (or dKP equation): 

(AUt - l2UUx)x = 3C/yy . 

This too is integrable, but the description does not use pseudo-differential operators but 
a more geometrical description in terms of vector fields and differential forms. 

The central idea of this paper is the development of a deformed differential calculus 
based on an associative T*r-product and its application to the theory of integrable systems. 
One will obtains an elegant description of the KP hierarchy in terms of vector fields and 
differential forms rather than the more usual pseudo-differential operator formalism. The 
advantage of this approach will be two-fold: firstly, one retains a geometric description, 
so, for example, one can go to a dual description in terms of differential forms; secondly 
in the limit in which the deformation disappears one recovers the dKP equation directly 
without the need of the somewhat singular limit outlined above. 

This deformed calculus will be derived in section 2 and used in section 3 where various 
examples of multidimensional integrable systems (an integrable deformation of the (anti)- 
self-dual vacuum equations, the KP hierarchy and the Toda hierarchy) will be studied. It 
will turn out the all these systems may be written in terms of a 2-form 17 satisfying the 
equations 



= 0, 
OAO = 

in analogy to Gindikin's bundle of forms. In section 4 the geometry of the KP hierarchy 
will be examined in more detail. This work raises a number of further questions, some of 
which are outlined in section 5. 



2. Deformed Differential Geometry 

A Poisson manifold M. is endowed with a bilinear skew-symmetric Poisson bracket defined, 
for u,ve C°°(A4), by 

and with the additional property that 

{{/, g}pB, h}pB + {{g, h}pB, f}pB + {{h, f}pB, g}pB = , 

this being known as the Jacobi identity. Further, it will be assumed that is a symplectic 
manifold, that is a Poisson manifold for which the matrix o;*^ is of maximal rank. It follows 
that the dimension of A4 must be even, so 



dim {M) = 2N 



for some integer . It will be assumed that uo^^ is constant and that a basis has been 
chosen in which 



ij _ I In 
[ -In 



This structure may be used to define a deformation of the above Poisson bracket. For 
u,v E C°°{J^) one defines a new product 

U'kV = exp f — Co'*-'-— -^:r-)tt(x)v(x) 

or, on expanding the exponential 

Ui,V=} tUj'^^\ . .UJ''^'— TT—TT-- TT—- (2) 

2^*3! . . . dx^^ 9x^1 . . . 9x> 

With this the deformed, or Moyal bracket, is defined as [Mo] 

{u,v} = . (3) 

K 

Lemma 1 

For constants c and functions u,v E C°°{Ai) , 

(a) lirrif^^Q U'kv = uv, 

(b) c-ku = cu , 

(c) -k is associative , 

(d) lim^^o {u, v} = {u, v}pb , 

(e) {u, v} is bilinear, skew-symmetric and satisfies the Jacobi identity. 
Proof 

Straightforward from definitions (P and (0). Note that the Jacobi identity follows 
from the associativity of the ^-product. 

□ 

The original motivation for the introduction of such a bracket came from a description 
of quantum mechanics using phase space variables [Mo]. Here k is replaced by —ih , and 
h —>■ corresponds to taking the classical limit. 

It is will necessary to introduce another product. For u,v E C°°{Ai) define 

^ 22^ (2s + 1)! ■ ■ ■ dx^^ . . . dx^^^ dx^ . . . dx^s ' ^ > 



Lemma 2 

For constants c and function u,v E C°°{Ai) , 



(a) u o V = V o u , 

(b) CO u = cu 

(c) lirrif^^o uo v = uv , 

(d) o is not associative, 

(e) 2d{n uov)/dK = u-kv + v-ku. 

Proof 

Again, these results follow from definitions and (^). 

□ 

With this o-product a deformed, or quantum, differential calculus will be constructed. 
A similar calculus has recently been constructed by Fedosov [Fe], this using the T*r-product 
rather than the o-product. The reason for the introduction of the new product will 
become apparent later (the motivation coming from the application of this calculus to 
multidimensional integrable systems) and rests on the following result: 

Proposition 3 

For u,v e C°^{M), 

du ^ dv _ 

Proof 

This follows from the definitions (§) and @: 

n — \ /.riJi ?,r2sj2s,,,'i 



00 - — o - — = > — — ttW . . . 00 ■'"■'■'''UJ 



dx'' dx' 22^(2s + 1)! ' ' ' dx^^ . . . dx^^-^dx'' dx^^ . . . dx^^^'dx' 



K ^ 22^+1 (2s + 1)! 9a;*i . . . ^x^^^'+i dx^^ . . . dx^^''+^ 



/U-k V — V -k U 
^ K 

{u,v} . 



In the limit /t — > this reduces to the definition (|T|). 



□ 



In the simplest case (A^ = 1) these formulae may be easily rewritten using the explicit 
form of 00^^ , coordinates x^ = x ^x^ = y : 



= E^E(-1)M J ]d:~'dludidpv, (5) 

oo 2s 2s / Oo \ 

""^ = g 2^.«(2. + 1)! g'^^' ( / J ■ (6) 

for functions u{x, y) , f (x, y) G C°°(A1) . 
Example 1 

Let = , the 2-torus. Functions on may be expanded in terms of basis 
functions 

Cm = exp i{mix + m2y) . 
With these one obtains from (0)-(0): 

Cm^en = exp(^m X n) Cm+n , 
{Cm, Cn} = Sinh^(m X n) Cm+n, 



sinh x n) 

m+n ) 



Cm o e„ = 2^^^ — 7^-^ e 



Km X n 



where m x n = m2ni — min2 . 



Given such a symplectic manifold and product one may define tangent and cotangent 
bundles TM. and T*Ai , the inner product between basis elements |^ and dx^ being given 
by 




The first difference is the formula for the inner product between general elements X G 
TpM and u G T*M , 



{X,u;) = {X^^,u;jdx^), 
= X'oui, 

i.e. the multiplication being done with the o-product. Similarly, given a vector field X 
and function / one defines 



A7 = X-o§, (8) 

again using the o-product. The general procedure should already be apparent: the only 
change to the standard, or undeformed, theory is when objects are combined, this being 
done with the o-product. Thus in the k ^ limit the standard theory is recovered. One 
may extend this new calculus to general tensor fields. However of more interest is the 
extension to an exterior differential calculus. 

The (i-operator, which maps r-forms to (r+l)-forms is defined as normal. For example, 
given a 0-form (i.e. a function) the 1-form df is defined by the relation 

{X, df) = Xf 
for all vector fields X . From this follows the formula 

df = ^dx' . 

The wedge product combines forms and so this will be done using the o-product. Explic- 
itly, if 



A = Ai^^i^dx'^ A ... A dx'" , 
B = Bj^,„j^dx^' A ... A dx^" 



then 



A A B = o Bj^_j^]dx'^ A ... A dx'" A dx^^ A ... A dx^" 



Example 2 

Suppose dim M — 2 (i.e. N — 1). Then for functions f{x, y) , g{x, y) e C°°(A<) , 

df = fxdx + fydy , 
dg = g^^dx + gydy , 

and hence 



df Adg = {f:c °9y- fy° 9x) dxAdy, 
= {/, g} dxAdy. 

Note that this uses the symmetry property of the o-product. Care must be taken in higher 
dimensions since, as the o-product is not associative, A A (B A C) 7^ (A A B) A C for 
arbitrary forms A , B and C . 



Having defined an exterior differential calculus one may define another intrinsic differ- 
ential object, namely a Lie derivative Lx corresponding to some vector field X e TM. . 
On functions 



and on vector fields 



This will also be written as CxY = [X^ Y] , and called the commutator of two vector fields. 
Using the symmetry of the o-product it follows that the commutator is antisymmetric. 
One may extend the definition to more general objects such as tensor fields in such a way 
that the theory is consistent. For example, for any p-form uj and vector field X , 

d{Cxoj) = Cx{duj) . 

Normally one has the relations 



[X,Y]f-X{Yf) + Y{Xf) = 0, 
[[X,Y],Z] + [[Y,ZIX] + [[Z,Y],X] = 0. 

The proof of these results uses the associativity of normal (i.e. undeformed) multiplication, 
and so do not hold for the deformed definitions based on the non-associative o-product. 
However for an important class of vector fields these results do hold. Given a function 
/ G C°°{Ai) (the Hamiltonian) the corresponding Hamiltonian vector field Xj is defined 

by 

Xf = uj-^^. (9) 

Strictly speaking these are local Hamiltonian vector fields, Hamiltonian vector fields hav- 
ing to be defined globally on M. . 

Lemma 4 

For functions f,g,hE C°°{A4) and corresponding Hamiltonian fields Xf, Xg and Xh , 

(a) X;/i = {/,/i}, 

(b) [Xf.Xg] = X{f^g} , 

(C) [Xf,Xg]h = Xf{Xgh)-Xg{Xfh), 

(d) [[Xf,Xg],X,] + [[Xg,X,],Xf] + [[X,,Xf],Xg] =0. 

Proof 



(a) This follows from definitions (|^), (|^) and Proposition 3. 
(b) 



— u 




■ dg d 




U,9} ■ 




[Xf,X,]h 



X{f,g}h, 

{{f,9}h}, 

{/, {9,h}}-{g,{f,h}}, 
Xf{X,h)-X,{Xfh). 



This uses result (b), the Jacobi identity for the Moyal bracket and the antisymmetry of 
the Moyal bracket (Lemma 1). 

(d) This follows from the Jacobi identity for the Moyal bracket (Lemma 1). 



Example 3 

For 2-dimensional manifolds A^^ these Hamiltonian vector fields generate the Lie alge- 
bra of area preserving diffeomorphisms of the manifold where the area element is dx A dy 
and the composition of two Hamiltonian vector fields is defined to be the Lie bracket of 
these fields. Explicitly, the field Xf generates the infinitesimal transformation x — > x—efy , 
y —>■ y + efx . This Lie algebra will be denoted by sdif f^iM^) ■ 



The differential objects constructed have been intrinsic to the manifold. One should 
be able to introduce a connection on A4 and define a covariant differentiation and hence 
curvature. A similar programme has been carried out using the ^-product by Vasiliev [V] . 
However for the application of this calculus to the theory of integrable systems such a 
structure will not be required. 

This ★-product is essentially unique. For a product 



□ 



oo 



f*9^ fg + Y.i^^'Qrif, g) 



r=l 



(where the Qr are bilinear differential operators) to be associative places strong restrictions 
on the type of higher-order terms that may be added. Further, the requirement that the 
bracket defined by {/, g}' = {f g — g -i< f) / n should reduce to the Poisson bracket in the 



«: — >■ limit introduces further restrictions and from these considerations follow various 
results on the uniqueness of the Moyal bracket [A,BFFLS,F1]. However these uniqueness 
results only state that any such deformations are equivalent to the Moyal bracket; there 
are apparently different structures which, after various changes of variable, become the 
Moyal bracket. For example one may define the following associative ^ir-product (in the 
N = 1 case) 

This too defines a bracket 

K 

which reduces to the standard Poisson bracket in the k — limit. This new bracket will to 
called the Kupershmidt-Manin bracket [K,Ma]. As above, one may define a corresponding 
o-product 

oo s s 

f-9 = Y: r— W E d'-'^d^f d^-^d^g , (12) 

s=0 l'^ + -"-J- m=0 



SO that 

du dv 



and hence an equivalent deformed differential geometry based on these new structures. 
The form of the ^-product is somewhat simplier then that given by (|), though the 
dependence on the sympletic structure of Ai is less transparent. The importance of these 
new products comes from their relationship to the algebra of pseudo-differential operators. 

A pseudo-differential operator P is an operator of the form 

finite 

p= y: «.(^)^' 

j=-oo 

where the multiplication of two such operators uses the generalised Leibnitz rule 



oo 



m(m - 1) ... (m - A; - 1) 



The set of such operators will be denoted by V . The symbol of a pseudo-differential 
operator is a function of two variables defined by 

finite finite 

sym{ ^ aj{x)d^) = J2 CLj{x)y^ . 

j=—OD j=—ao 

It has the important property that for all P ,Q E V 

sym{PQ) = sym{P) sym{Q) 



where *k=i denotes the ^-product (p!0| ) evaluated at k = 1 . It follows from this that 



sym{[P,Q]) = {sym{P),sym{Q)}'^^-^ 

where [P, Q] = PQ — QP . Thus one may replace pseudo-differential operators and its cor- 
responding algebra by functions of two- variables where the composition of two functions 
is done with the Kupershmidt-Manin bracket evaluated at k = 1 . More details of these 
algebraic properties may be found in [F-FMR] . Using the ideas developed above one may 
give these pseudo-differential operators a geometrical interpretation. 

Theorem 

Let Ti be the space of Hamiltonian vector fields (where = 1) whose Hamiltonians 
have Laurent expansions in the variable y , 

finite 

n = {Xj :/= Y: a,ix)y^}. 

j=-oo 

Then the map 

L : P/{constants} Ti 

given by 

l{P) = Xsym{P) 

is an isomorphism. Moreover 

'•([-P)<5]) = [Xsym{P)-,Xsym{Q)\ 
= ^{sym{P),sym{Q)y^^^ 

where the Lie bracket of vector fields is evaluated using the product o^^^i given by (|1^) eval- 
uated at K = 1 . 

Proof 

Straightforward. Given a Hamiltonian vector field one can construct the corresponding 
Hamiltonian (up to a constant) and hence a pseudo-differential operator whose symbol is 
the Hamiltonian. Conversely, l{P + c) = l{P) for all P eV . The last part of the theorem 
follows from the properties of the symbol map. 

□ 

The set of Hamiltonian vector fields clearly forms a Lie algebra under the composition 
defined by the Lie bracket. One may define the adjoint representation as follows. For 
functions f ,g,F e C°^{M) define 



ad{f)g = {f,g}, 
Ad{F)g = Fi<gi<F-^, 



the connection between the two coming from the deformed exponential 

oo 

exp« f = l + J2—f-k...-kf . 



n=0 

n— terms 



So, if F = exp^ / , 

Ad{F) g ^ exp^{ad{f)) g 
(this uses the Baker-Campbell-Hausdorff formula). On vector fields, 

ad{Xf) Xg = X{f^gy , 

Ad{F)Xg = Xp^g^F-l. 

These will be used in section 4 to describe the dressing properties of the KP hierarchy. 
The residue of a pseudo-differential operator P — J2 dnd^ is defined by 

res{P) = a_i . 

It follows that 



res{P) — res{sym{P)) , 



where the residue of the function sym{P) is the normal residue, regarding sym{P) as 
a complex function of y . The residue has many uses, in particular in the study of the 
Hamiltonian properties of integrable systems. 



3. Applications to Integrable Systems 

In this section a number of multidimensional integrable systems will be studied using the 
geometric structures developed in the last section. It will be shown that these systems 
may all be written in terms of a 2-form f2 which satisfies the equations 



rfri = 0, 

OAO = 0. 

These equations encode the integrability conditions for these systems in an elegant geo- 
metric manner. 

Let A1 be a sympletic manifold with some associated ^-product. In applications one 
will require an extended manifold 



M = MQ)r 



where T consists of an extra set of coordinates (for example the 'times' in a hierarchy of 
evolution equations). The manifold Ai may be thought of as a phase space and in the 
applications considered here this will be two-dimensional. The T*r-product extends to a 
product on A4 by 

d d 

«(x, t) -kv{^, t) = exp (/€ J^ ——)u{^, t)t;(i, t) |^^. 

(where x = {x*} are coordinates on M. and t are coordinates on T), that is, the depen- 
dence on the coordinates on T is ignored. The differential calculus outlinded in section 2 
similiarly extends to the manifold M. . One difference is that Xj will refer to a Hamiltonian 
vector field on M. whose Hamiltonian function depends on the coordinates on M. , 

_ dfi^d 

Af — UJ — 7- — r , 

i.e. a time dependent Hamiltonian vector field on Ai , where the 'times' are the coordi- 
nates on T . 



3.1 The Anti-Self-Dual Vacuum Equations 

The anti-self-dual vacuum equations govern the behaviour of complex 4-metrics of signa- 
ture (-I-, -|-, -|-, +) whose Ricci tensor is zero and whose Weyl tensor is anti-self-dual. Since 
these curvature conditions are invariant under changes of coordinates there are many ways 
to write these equations. On particular form of the equations uses the fact that such met- 
rics are automatically Kahler and so may be written in terms of a single scalar function 
Q , the Kahler potential. The curvature conditions then give the equation governing the 
potential (known as Plebanski's 1^* Heavenly equation [PI]): 



dxdx dydy dxdy dydx 
The corresponding anti-self-dual Ricci-flat metric is 



(13) 



g{^) = -^—^dx'dx\ x' = x,y, x^=x,y. (14) 

This equation can, in principle, be solved using a Penrose transform - the original non- 
linear graviton construction [Pen]. Although not realised at the time, the existence of 
such a transform makes (|l^) into a completely integrable system, an important, and rare, 
example of a multidimensional integrable system. As such it has all the properties one 
would expect of such a system, an infinite number of conservation laws [S93] and an 
associated hierarchy [S95b], for example. A Lax pair for this equation was derived by 
[NPT] and later interpreted by Park [Pa] as that for a 2-dimensional topological chiral 
model with gauge potentials in the infinite dimensional Lie algebra sdif f{M.'^) for some 
two dimensional manifold M. . 



The equation may be written as 
where the Poisson bracket is defined by 

{/; fl'jpB = fxdy — fy dx ■ 

The equation that will be studied in this section is an integrable deformation of this 
equation, where the Poisson bracket has been replaced by the Moyal bracket (0) 

{n^,ny] = i. (15) 

The space M. will have coordinates {x, y} (with composition using the products (H) and 
(|^)) and the space T will be taken to be (or possibly C^) with coordinates {x, y} . This 
deformed system ([TsD will be studied using the deformed calculus developed in section 
2. It will turn out the it shares many of the features and properties of the undeformed 
system (plsl) . 

Let U and V be the following vector fields on TM. 



d_ 



V = \— + Xg, 



where A G CP^ is a constant known as the spectral parameter. The system of equations 
for the function ip e C°°{M) 



m) = 0, 

V(V') = 



(or, equivalently. 



X^y + {g,^} = 0) 
is overdetermined unless the integrability condition 



[U,V]=0 

holds. Here [, ] is the Lie bracket of vector fields. If these equations are satisfied, then 
one has two independent solutions L and M for ip which satisfy the equation {L, M} = 1 . 
Note that the above equations may be written in the following ways: 



U{L)=U{M) = 0, 
V(L)=V(M) = 0, 



or 



(W, dL) = {U, dM) = , 
(V, dL) = (V, dM) = . 

The above integrability conditions are satisfied if tfie functions / and g satisfy the equa- 
tions 

fy-gx = , 
{f,9} = 1. 

The first equation implies the existence of a scalar function Q such that f = , g = 
and with these the second equation becomes the deformed Plebanski equation (piSl) . This 
shows that this may be interpreted as a 2-dimensional chiral model with gauge potentials 
in the Lie algebra sdif fi^{A4'^) . 

In [S92] the vector fields U and V were interpretated as operators: 

oo 2s+l 2s+l / Oo I 1 \ 

^ = + E ^.iT^^j^ E (-1)^- ( ) ^^'-m m^'-' ^ 

oo 2s+l 2s+l / 9„ I 1 \ 

V = Aa, + g ^...fp^ ^ ^, S (-1)' ( + J ar-'aggalar'-^ . 

The geometrical approach used here is much simplier, and the manipulations using the 
o-product which lead to equation (^) are almost transparent, deviating very little from 
the undeformed calculation which leads to equation ([13D (to achieve such a result was 
one of the original motivations in the development of the deformed calculus). Another 
advantage of the geometrical over the operator based approach is that one can go over to 
a dual description in terms of differential forms on T*(A^) . 

Let Vt be the 2-form 

rt = dx /\dy + XiVL^xdx /\dx + VL^ydx Ady + flyxdy Adx + VLyydy A dy) + X^dx A dy . 
This clearly satisfies the condition dVt = , and in addition 

ri A ri = X'^{flxS: ° ^yy — ^xy ° ^yx — i)dx A dx A dy A dy , 
= X^{{Vtx,VLy} -l)dx Adx Ady Ady , 
= 

by virtue of (p!5|). Thus the Lax pair, and hence the integrability of this deformed system 
is encoded into the equations 



dri = 0, 

nAn = 0. 



Further properties of this system have been found. In [S92] a perturbative solution in 
powers of k was constructed. On writing 

oo 
n=0 

one obtains Plebasnki's equation ([I3D for and an infinite number of hnear equations 
for the f2„ , n > of the form 

^Uo^n = Sni^O,--- , ^n-l] , 72 = 1, 2,. ..,00. 

The operator Dq is the wave operator on the spacetime given by the metric g{^) given 
by equation (|14D and the function Sn is some known function of its arguments. A similar 
procedure may be applied to a Moyal-algebraic deformation of Plebanski's 2"^^^ heavenly 
equation [PPRT]. In [C] a symmetry reduction of this system was studied and in [T94a] 
the dressing transform (using a Riemann-Hilbert factorisation in the corresponding Moyal 
loop group) was constructed. As mentioned earlier, one may study the conservation 
laws, symmetries and hierarchies associated with Plebanski's equation and these results 
still holds, under the replacement of the Poisson bracket by the Moyal bracket, for the 
deformed system ([T5|). 

A slightly more general framework may be achieved by observing that the vector fields 
dx ,dy ,Xf and Xg which make up the vector fields U and V all preserve the volume form 
Lj = dx A dy A dx A dy on , as m the construction of self-dual metrics [MN]. 



3.2 The KP hierarchy 



The KP hierarchy is defined as follows. Let C be the pseudo-differential operator 

oo 

/: = a+^M„(x,t)a-" 



n=l 



where t = {ti , t2 , • • •} • The evolution of the fields Un{x, t) with respect to the times t is 
given by the Lax equations 

f = |B„.£| (16) 

where 

i3„ = [£"]+, n = l,2,...,oo 

and [C]+ denotes the projection onto the purely differential part of the pseudo-differential 
operator O . Similarly, denotes the projection onto purely negative powers, so O = 

o+ + o^. 

Let the coordinates on M. be {x. A} and times t be coordinates on T . The -k and o 
products on M. will be given by equations ([T0| ) and (|12D evaluated at k = 1, and the ' will 
be dropped on the corresponding bracket for notional convenience. Taking the symbols 
of the above operators gives the following functions on M. : 



L = sym{C) , 

oo 

= A+ ^M„(a;,t)A-", 

n=l 

Bn = sym{Bn) ■ 

Note that 

sym{C) = sym{C) *k=i . . . -k^^i sym{C) , 

(which does not equal sym{C)"-) which Kupershmidt [K90] denotes by L*" (see also [F- 
FMR]). The Lax equation ([l6| ) becomes the vector field equation 

L„(L) = 0, n=l,2,...,oo. (17) 

Here L„ G TA4 is the vector field 

and the operations are performed with k = 1 . At this point this condition will be dropped, 
thus obtaining a ^-dependent KP hierarchy. The Lax function L remains unchanged, but 
the Bn acquire k dependence, since they are now defined by the formula 

(and so reduce to the previous definition if k = 1), where + denotes the projection onto 
non-negative powers of A, . This has the advantage that one may recover the dispersionless 
KP hierarchy in the /t ^ limit without the need for rescaling variables. 

Example 4 

The first few equations in the ^-dependent KP hierarchy ( |17D are [K90] 



Bi = X, 
B2 = X^ + 2u2 
B3 = + SXui + 3u2 + 3kui . 



which leads to the evolution equations 



M2,t2 = 2M3,a; + 2MiMi,^ + KU2,xx , 

Mi_t3 = 3M3,a; + 6MiMi_a; + 'iKU2^xx + I^^Ui^xxx 

(the ti-fiows are trivial). These show the ^-dependent terms. On eliminating U2 and M3 
one obtains the KP equation 



One may obtain an equivalent KP hierarchy by using the Moyal ^ir-product and bracket 
rather than the Kupershmidt-Manin T<r-product and bracket [K,S95a]. One obtains the 
functions 



Bi = A, 

= + 3Ami + 3u2 
(^-dependent terms only appear in the -B„ for n > 3) and evolution equations 



Ul,t2 = 2^2,0; , 

Ml,t3 = GUiUi^x + 3U3,X + -^Ui^xxx ■ 

This system also leads to the KP equation on eliminating U2 and M3 , and on redefining the 
fields it is easy to see that these two systems are equivalent. Note that in both cases the 
limit K — i> one obtains the dispersionless KP equation directly without further rescaling 
of the variables. 



An alternative form of the KP hierarchy is based on the zero-curvature conditions 
(which follow from the Lax equation (|16|)) 



dt dt 
or equivalently 

dBn OB 

m 



+ 5„} = . (1^ 



Note that this is the condition of the vector fields L„ to commute, [L^, L„] = for all 
m,'n, = l,2,...,cxD. 

These zero- curvature relations may be encoded in a 2-form O defined by 

00 

Q, = d\ Adx + ^ dBn A dtn ■ 

n=2 

This form satisfies the following equations 



dfl = 0, 
riAfi = 0. 



The first is obvious. One has to be careful in evaluating the second equation (see Example 
2 in section 2), but one obtains 



m,n=2 ^^'n% ULn 



dX A dx A dtm A dty 



and hence A = if and only if the zero-curvature relations (|1^) hold. The geometry 
of the KP hierarchy will be discussed further in section 4. 



3.3 The Toda hierarchy 

The definition of this hierarchy is very similar to the definition of the KP hierarchy. The 
Lax operator is 



n=0 



nd 



(note the range of summation) and the evolutions of the fields are given by the Lax 
equation 

where 

S„ = [/:"]+, n = l,2,...,oo 

and [0]+ denotes the projection onto positive powers of of the pseudo-differential 
operator O . The operator acts as a shift operator, 

e'^^/(x) = /(x + n). 

The geometric description of the hierarchy is obtained in the same way as above. 
Taking symbols of the operators give 



L = sym{C) , 

oo 

Bn = sym{B„) 
and the above Lax equation becomes 



{Bn,L}. (19) 



dL_ 

dtn 

Once again the condition k = 1 will be dropped, so composition will be done using the 
Kupershmidt-Manin -k- and o-products (|10|) and (|12D, so now the are defined by the 
equation i?„ = where + denotes the projection onto non-negative powers of , as 

in the K-dependent KP hierarchy. In the limit k ^ one obtains the dispersionless Toda 



hierarchy. One difference between the hierarchy and the KP hierarchy is that the evolution 
equations for the fields contains an infinite number of /t-dependent terms. However these 
may be recombined in terms of shift operators, as the following example will show. 



Example 5 

One possible truncation of this hierarchy is to set m„ = for n > 2 , so that 

L = + uq + uie~^ , 

The evolution equations for the fields uq and ui are given by 

dL 



dt 



where the bracket is the Kupershmidt-Manin bracket (^) and, for greater generality, the 
K = 1 condition has been dropped. 
This gives the equations 



ls=0 



Ul{x) 



Ui{x + K) — Ui{x) 
K 



K 



s=0 



Uq{x) , 



Ui{x)[uq{x) - Uq{x - k)] 



and on eliminating uq on obtains the Toda lattice equation 

uAx + k) — 2ui{x) + uAx — k) 



(log Ml (x) )tt 



Instead of using the Kupershmidt-Manin bracket one could use the Moyal bracket. 
This gives the slightly different equations 



ui,t{x) 



K 



K 



y — 92^+^ 



ui{x) 



Ui{x + k/2) — ui{x — k/2) 



2ui{x) 



K 



\y I 



2s+l 



s=0 



22^+i(2s + 1)! 



i2s+l 



uo{x) 



Ui 



uo{x + k/2) — uq{x — k/2) 



On eliminating uq one again recovers the Toda lattice equation. Note that with the 
Kupershmidt-Manin bracket one obtains a forward/backward difference operator while 
the Moyal bracket gives a central difference operator. In both acts as the lattice 

spacing and as k — one recovers the dispersionless Toda equations (since both these 
brackets are deformations of the Poisson bracket) 

Uo,t = Ui^x , 

Some of the properties of this system and its hierarchy may be found in [K85,FS]. 



As with the KP hierarchy, the Lax equation ([19D is equivalent to a set of zero- curvature 
relations for the and these may be encoded into a 2-form f2 which satisfies the equations 

dn = 0, 
flAfl = 

in exactly the same way as was done for the KP hierarchy. 

4. The Geometry of the KP hierarchy 

The main result of this section is to show how a solution of the KP hierarchy may be 
associated to a solution of a Riemann-Hilbert problem in the Lie group SDif fi^{M.'^) (the 
Lie group corresponding to the Lie algebra sdif fi^{A4'^) . Explicitly, given a map 

( x] ( f{x,k) \ 
[k) [g{x,k)) 

with {f,g} = 1 then this map factors, so there exists a map 

[q) [9{P,Q)) 

where the right hand side is analytic in k (the notation S- will be used to denote the 
part of the Laurent series S consisting of negative powers of k only). The results derived 
in section 2 enable existing results on the KP hierarchy to be lifted whilst furnishing 
them with a geometrical interpretation (the definitions of the manifold Ai and -k- and 
o-products will be the same as in section 3.2). In this section k, = 1 and the K-symbol on 
the exponential exp^ will be dropped. The main results of this section are due to Takasaki 
and Takebe [TT] . A more careful analysis is needed for n ^ 1 . 

More fundamental than the Lax operator C is the operator W defined by 

oo 

n=l 



with which the Lax operator is defined as 

The evolution of the fields w„ are given by the equation 

i3„W - Wa" 



from which follows the Lax equation (piq) . On taking the symbols of the operators one 
obtains 



W = sym{yV) , 

oo 
n=l 

L = sym{C) , 
= Ad{W)k. 

The Orlov operator Ai is defined by [GO] 

oo 
n=l 



or, equivalent ly, by 



M = sym{M) , 

= Ad{W e^^[t{k)\)x 



where t(A;) = E^=itn/^" • 
Lemma 

The pair (L, M) satisfy the equations 

{L,M} = 1. 

Conversely, given such a pair then there exist a unique dressing function W so that 
L = Ad{W)k and M = Ad{W exp[t{k)] )x . 

□ 



Such a (L, M) pair will be said to satisfy the KP hierarchy. The first two of these 
equations may be re-written as vector field equations: 



{dt„-XBjL = 0, 
{dt„-XsjM = 0, 

or alternatively, using the inner product (, ) between vector fields and forms, as 



{dt^-XB„,dL) = 0, 
{dt„-XB„,dM) = 

for n = 1 , 2 , . . . , oo . These functions (L, M) will play the analogous roles to the coordi- 
nates on the twistor surfaces in the nonlinear gravition construction. 

The next theorems show how such a pair are related to a Riemann-Hilbert factorization 
problem. The first show how a solution to the Riemann-Hilbert problem defines a solution 
to the KP hierarchy and the second shows the converse. 

Theorem 



Suppose on has functions 



L = Ad{W)k, 
M = Ad{W ex.p[t(k)])x 

(with {L, M} = 1 ). Then for any pair of functions f{x, k) , g{x, k) (with Laurent series) 
if 



{f,g} = 1, 

/(M,L)_ = 0, 
g{M,L)_ = 0, 

then the pair (L, M) satisfies the KP hierarchy. 

□ 

Theorem 

If the pair (L, M) satisfy the KP hierarchy then there exist functions / , g such that 



{f,9} = 1, 
/(M,L)_ = 0, 
g{M,L)_ = 0. 



The proofs are basically identical to the proofs in [TT], the only difference being that 
here they are reformulated in terms of the deformed geometric structures rather than in 



terms of pseudo-differential operators. One may also prove an uniqueness results for the 
solution (L, M), at least in the neighbourhood of the trivial solution {k, x) . 

One outstanding problem is how to relate the 2-form O to the pair (L, M) . In the 
dispersionless limit one has (for all the systems discussed in section 3) a relation 

n^dLAdM. 

However, the proof of this result relies on the associative property of normal multiplication 
which nolonger holds for the deformed o-product. It may be that this result still holds, for 
example, the obstruction might vanish due to the relation {L, M} = 1 . This problem, of 
how to understand the direct relation between the pair (L, M) and fl is currently under 
investigation. 



5. Comments 

In summary, the three classes of integrable hierarchy discussed in section 3 may all be 
formulated in terms of vector fields which preserve a volume form 

oo 

u! = dx f\dy f\ l\ dtn 

n=l 

(or uj = dx Ady Adx Ady in the deformed anti-self-dual vacuum equations) on 7W , together 
with function (L, M) which satisfy the equations 

V^{L) = 0, 

Vi{M) = 0, 

or, in terms of the inner product between TM. and T*M. , 

(V„rfL) = 0, 
{Vi,dM) = 0, 

with the functions L and M being independent: {L, M} = 1 . A dual description also 
exists for all these systems in terms of a two form Q on Ai satisfying the equations 

dn = 0, 
flAfl = 0. 

The precise relationship between these two dual descriptions requires further investigation. 
In all cases the solutions are encoded in a Riemann-Hilbert problem in the corresponding 
loop group. 



This work raises a number of further question. For example, it should be straight- 
forward to examine the sysmmetries of these integrable systems using thses methods, 
with the Hamiltonian vector fields playing the role of the symmetry generators (see, for 
example [T94b]). This would provide a geometrical description of various Wq^ and Wkp 
algebras. One use of such symmetries is in the construction of a constrained KP-hierarchy. 
One example of this contains the KdV hierarchy. However, the KdV hierarchy has been 
shown to be a reduction of the self-dual Yang- Mills equations (and its generahsations) . 
Thus there are two way of looking at the KdV equation: one based on the Yang-Mills 
self-duality equation and one based on the deformed differential geometry constructed in 
section 2. Precisely how these two seemingly different constructions are related deserves 
further study. In connection with this is how to understand the non-local nature of the 
Riemann-Hilbert problem for the KP equation compared with the local one for the KdV 
equation [AC,M]. 

All of the examples of integrable systems studied in this paper have used Hamiltonian 
vector fields in their construction. Are there any systems which use more general, non- 
Hamiltonian, vector fields? The property of commuting flows for these hierarchies can be 
traced back to the Jacobi identity for Hamiltonian vector fields, so any hierarchy based 
on non-Hamiltonian vector fields might lose this property. 

One possible use of this deformed calculus would be to develop a theory of deformed (or 
quantum) twistor spaces (which would encode the Riemann-Hilbert problems in the Moyal 
loop group) more axiomatically. One obvious place to start is to deform the symplctic 
structure on the fibres of the nonlinear graviton's twistor space. An observation that 
might be of use is that ^^r-product do exist on the complex manifold CP^ and other 
complex coset spaces. This suggest that one should develop a deformation theory (in the 
sense of Kodaira) of such spaces. Such idea, however, are outside the scope of this paper. 
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Notes added in proof 

Since this paper was written a number of other papers have appeared. In [PP] and [G- 
CPP] the Moyal deformation of self-dual gravity has been studied using a chiral model 
approach and in [S96] it was shown that the Toda lattice is a reduction of this Moyal 
deformed self-dual gravity, a result analogous to the reduction from the standard, unde- 
formed, self-dual gravity equations to the Boyer-Finley equation. Other notable papers 
are [Ke], [KeS],[DM-H] and [R], which develop various connections between discrete sys- 
tems, geometry and associative ★-products. 
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